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. Abstract 

> , 

f" — , We study the possibility of measuring nonstandard CP violation effects through Dalitz plot 

analysis in B ± — > D^JJ tt decays. The accuracy in the extraction of CP violating phases is 
analyzed by performing a Monte Carlo simulation of the decays, and the magnitude of possible 
new physics effects is discussed. It is found that this represents a hopeful scenario for the search 
C*~) ■ of new physics. 



CO 

O 



O 
43 



Oh" PACS numbers: 11.30.Er, 12.60.-i, 13.25.Hw 



• I Introduction 



The origin of CP violation in nature is presently one of the most important open questions in 
particle physics. Indeed, the main goal of the experiments devoted to the study of B meson decays 
is either to confirm the picture offered by the Standard Model (SM) or to provide evidences of CP 
violation mechanisms originated from new physics. Among the various CP-odd observables in B 
physics, attention is mostly concentrated in the "gold-plated" channel B — * J/\t K$- According 
to the SM picture, from the analysis of a time-dependent CP asymmetry observed in these decays 
it is possible to get a "clean" measurement of sin 2/3, where (3 is one of the angles of the so-called 
unitarity triangle pQ. Recent measurements by BELLE and BaBar Collaborations, together with 
previous results from Aleph, Opal and CDF, lead to the (averaged) value sin 2(3 = 0.734 ± 0.054, 
which is in good agreement with the constraints imposed by other measured CP-conserving and 
CP-violating observables [2]. 

In fact, the common belief is that the SM is nothing but an effective manifestation of some 
underlying fundamental theory. In this way, all tests of the standard mechanism of CP violation, 
as well as the exploration of signatures of nonstandard physics, become relevant. One important 



characteristic of the SM is that it includes only one source of CP violation, namely a complex 
phase in the quark mixing matrix Vqkm- In general, since overall phases of transition amplitudes 
are not observable, one has to deal with interference effects in order to search for measurable CP- 
violating quantities. Within the SM, there are some specific processes in which the amplitude is 
either dominated by a single contribution, or in which several contributions are significant, all of 
them carrying the same weak phase. In these situations, weak SM phases are unobservable, and 
asymmetries between CP conjugated processes are expected to be vanishingly small. This offers an 
attractive window to search for evidences of new physics, and is the main motivation for this work. 

We show here that three body decays B + — > DfD°ir° and B~ — » D~D°ir° provide an interesting 
scenario to look for such effects. For these processes, the main contributions to the decay amplitude 
in the SM carry a common weak phase, therefore the measurement of relative CP-violating phases, 
leading to an asymmetry between B + and B~ decays, would represent a signal of new physics. We 
discuss here the possibility of performing these measurements by means of a Dalitz plot (DP) fit 
analysis. In general, three body decays of mesons proceed through intermediate resonant channels, 
and the DP fit analysis allows a direct experimental access to the amplitudes and phases of the main 
contributions |Hj • The usage of this technique for a clean extraction of CP-odd phases has already 
been proposed in the literature [H] in relation with other CP-violating observables, more precisely, 
to get clean measurements of the weak angle 7 within the SM. From the experimental point of 
view, the usage of charged B mesons has the advantage of avoiding flavor-tagging difficulties. In 

addition, the processes B ± -> DfD° 7T° appear to be statistically favored, in view of their relatively 
high branching ratios of about 1%. 

In order to evaluate the experimental perspectives, we perform a Monte Carlo simulation of the 
actual processes, applying the DP fit technique to evaluate the error in the extraction of possible CP- 
violating phases. Then we perform a rough theoretical analysis, discussing the expected magnitude 
of CP violation effects that could arise beyond the SM. According to our results, the considered 
channels offer a promising scenario to obtain a clear signature of new physics. In the worst case, 
the lack of evidences would allow to improve the present bounds on the parameters of the model 
under consideration. 

The paper is organized as follows. In Sect. II we describe the general framework, introducing 
the CP-violating observables. In Sect. Ill we detail the DP fit procedure and present the results 
of our simulations. Sect. IV is devoted to the theoretical discussion of new physics effects, while in 
Sect. V we summarize our main results. 

II CP-violating phases and Dalitz plot fit technique 

In this section we describe how the DP fit technique can be applied to disentangle possible effects 

of new physics in B^ — > D^D Q ir Q decays. In principle, these processes are expected to proceed 
through various intermediate resonances, as well as through a direct, nonresonant channel. The 
total branching ratio will result from the interference of all these contributions. The Dalitz plot 
maximum likelihood technique is a powerful tool to get a clean disentanglement of the relevant 
contributing channels, allowing to measure the ratios between the different partial amplitudes 
together with their relative phases. This can be used to perform a clean extraction of CP-violating 
phases, avoiding many theoretical uncertainties. 

Let us begin by summarizing the main steps of this procedure. More details on these ideas 
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can be found in Refs. [BJ. In general, for a given three body decay, in the DP fitting analysis of 
experimental data one defines a fitting function J-{rn\ 1 m'^), where m\ and m| are the usual DP 
phase space variables. In our case this function can be written as 

^B^D a DO^{m\,ml) = l^jaje 1 ^ Ajimj^l)] 2 , (1) 

where m\ = (p n o +p D o) 2 , m 2 , = (p^o +p Da ) 2 , Aj are definite functions corresponding to each partial 
channel, and aj and Oj are real parameters that emerge as outputs from the fit. The index j labels 
the intermediate resonant channels, as well as the nonresonant one. For the resonant channels, 
the main phase space dependence of the functions Aj is given by the Breit-Wigner (BW) shape 
characterizing the resonances, together with definite angular functions which depend on the spin of 
the corresponding resonant state (we come back to this issues in the next section). The nonresonant 
decay amplitude is assumed to be constant in most experimental analyses. This fitting technique 
has proven to be successful to describe e.g. three body decays of D mesons (7j- In those analyses 
the phases Oj have been extracted with combined statistical and systematic errors as small as a few 
degrees, in experiments with a few thousands reconstructed events. 

In general, the phases Oj can be written as the sum of a "strong" (CP-conserving) phase 5j and a 
"weak" (CP-violating) phase <pj. These cannot be measured separately by a single fit. Nevertheless, 
comparing the outputs from the CP-conjugated B + and B~ decay experiments one can extract 
both phases 5j and tpj simply from 

*i = \Wt + 6 J) ( 2 ) 

<Pj = \( e t~ e J)i ( 3 ) 

where 0j~ (Oj) stands for the phases measured from the B + (B~) decays. It is worth to notice that 
weak phases can be extracted even in the limit where strong phases 8j are vanishingly small — which 
is expected to be the case in many B decays, owing to the large b quark mass. This represents a 
remarkable advantage with respect to most proposals of measuring CP asymmetries in charged B 
decays. In general, in order to get a sizable asymmetry, one requires the presence of strong FSI 
phases, which introduce a significant theoretical uncertainty. In our case, however, strong phases 
are already supplied by the resonance widths in the BW functions 9 u , and no theoretical estimation 
of FSI phases is needed. Moreover, the latter can be independently obtained from the fit by means 
of Eq. Q. 

It is important to notice that for the fitting procedure to apply, it is necessary that the decay 
amplitude receives contributions from at least two intermediate channels carrying different CP 
violating phases. Indeed, an overall phase is physically meaningless, and the DP fit only allows the 
measurement of relative phases between different channels. 

Let us now analyze the case of the decays B ± -» DfD° 7T° in the framework of a theory including 
physics beyond the SM. For each intermediate channel, it is natural to assume that new physics 
occurs at a relatively high energy scale, therefore its effects can be decoupled from the resonance 
BW functions, the angular functions, and other possible form factors in Aj{m%, m|) arising from 
low-energy hadronic interactions. Accordingly, the fitting function T will be still of the form 
proposed in Eq. (|Tj). with the same functions Aj, and now the complex weights otj e j6,J will include 
the effects of new physics: 

±(exp) ie ±(exp) SM i(8 SM ±<p SM ) . NP i(S NP ±>p NP ) , ,s 

e 3 = aj e y J r i ' + a j e ^ ^ ; . (4) 
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Here the index (exp) refers to the experimentally measurable quantities (outputs of the DP fit), 
whereas SM and NP denote Standard Model and new physics contributions respectively. The ± 
signs correspond to decays of mesons. 

Within the SM, the short-distance effective Hamiltonian relevant for the decays — ► D^D ^ 
has been studied in detail ^01, including the effects of strong and electroweak penguin operators. 
The situation is similar as in the "gold-plated" channel B — * J/^/Ks, in the sense that the main 
contributions (both tree and penguin) to the effective operators carry the same weak phase. In 
this way, this phase is expected to factorize, being common to all (resonant and nonresonant) 
channels contributing to the decay. One could thus conventionally set ipj M = for all j, and no 
CP asymmetry should be observed between B + and B~ decay patterns in absence of new physics 1 . 
On the contrary, if new physics is present, the situation may be different. To simplify the analysis, 

let us assume that only two intermediate channels contribute, namely those mediated by resonances 

H 

D* v and Dg — say channels 1 and 2, respectively. This is a natural assumption, since in fact they 
are expected to largely dominate the decay (in any case, if other intermediate channels were shown 
to provide significant contributions, the procedure we describe here can still be followed on the 
same grounds). In this two-channel case, only the relative phases between both contributions 1 
and 2 can be measured. From Eq. (|4"|). one has 



d ±(exp) _ ±(exp) _ e ±(exy) = ^ 



a fM + a NP e i { 8r-5f M ±^ P ) 



+ 5f M - 6» M . (5) 

The theoretical framework can be simplified by introducing some natural assumptions. First, it 
is reasonable to think that the resonance hadronization and decay processes — which are governed 
by strong interactions in the nonperturbative regime — can be disentangled from the effects of 
new physics, the latter taking place at a high energy scale. In addition, it is usual to assume 
that strong FSI are the main source for strong phases 6j, since high energy contributions to CP- 
conserving phases arising from absorptive parts of QCD and electroweak loop diagrams are shown 
to be suppressed ^2j. In this way, for each resonant channel strong phases should factorize out, 
i.e. 5j = 5j M = 5j. On the other hand, in most scenarios of new physics — as well as in the SM 
itself — , CP-violating phases are essentially determined by the flavor content of the quarks entering 
the diagrams that dominate the b quark decay. If this is the case, since both resonant states 
D*°D S and D°D* have the same quark content, one expects that the new CP-violating phases 
obey ifi = (f2 P — l f NP '> remaining constant along the phase space. Once these assumptions have 
been taken into account, the measurable complex weights in Eq. © can be written as 

(exp) J9f^ _ ,SM , NP p ±i^h JS. 



a 3 e 3 



(a^+afe^ >) , j = 1, 2 . (6) 



We have dropped here the ± signs in aij , since the assumption 5^ p = 5j M implies a- 
( e xp) 2_ rpjjg expression for the relative phases Q^( ex P) in Eq. (JSJ) simplifies now to 



SM 



* ±( -'^ Q M ;>;^ )+*-*■ CD 



1 Within the SM, one expects in fact a tiny CP asymmetry in the decay rates. This has been analyzed in Ref. 
for the process B~ — > D°Dj , where the effect is found to be about 0.2%. 

2 In principle, this last relation could be experimentally checked through the comparison between the fits for B + 
and B~ decays, providing a consistency test for our assumptions on the strong phases. 
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As stated, we are interested in the difference between the relative phases for the CP-conjugated 
decays B + and B~ , which is an observable of CP violation. This is given by 

Finally, assuming that new physics contributions are small when compared to SM amplitudes, i.e., 
a^ p <C ctj, we end up with 

a? p a NP 



a 5M ft SM 



Notice that this quantity is independent of any CP-conserving phase. It only depends on the (in 

%I M and af 



principle, unknown) real amplitudes a^ M and af p , and on the (also unknown) CP-violating phase 



NP 

if 



As we have discussed above, A9^ exp ^ vanishes in the absence of new physics. This is a convenient 
situation for the search of clean effects of physics beyond the SM, provided that the factors in the 
r.h.s. of Eq. Q are large enough to allow a clear experimental signature. The perspectives in this 
sense are addressed in the next sections. 



Ill Experimental perspectives 

In this section we make an estimate of the precision that may be reached in the measurement of 

the phase difference A0( exp ^ in — > D^D ^ decays. This will indicate, according to the result 
in Eq. Q, the minimum size of new physics contributions to the decay amplitudes needed to yield 
a distinguishable experimental signal. 

One important reason for which the channels considered here deserve special attention is their 

relatively high statistics. Since the branching ratios for B^ — > D^D°ir are as large as ~ 1%, after 
a couple of years of full run of LHCb, and assuming a 20% reconstruction efficiency, one should end 
up with some 10 5 reconstructed events in each B + and B~ Dalitz plots. This is a large number, 
taking into account that DP fits performed for D meson decays with much less events have led 
to the measurement of relative phases with statistical errors of just a few degrees [Hj. However, 
the processes considered here are very different from those. Indeed, even if such a large number 
of events will certainly give a very precise measurement of the branching fractions for each partial 
channel — i.e., the quantities a^ exp ^ — , a precise measurement of phases requires not only large 
statistics but also a large interference region between the different intermediate channels. It is not 

obvious that this will be the case for B ± -> DfD°TT° decays, since the involved resonances are very 
narrow, their widths laying below 1 MeV |Sj. 

In order to evaluate the actual experimental feasibility of our proposal, we have carried out a 
Monte Carlo simulation of the decays. Our goal is to generate 10 5 events in the Dalitz plot, and 
then to perform a Dalitz plot fit analysis in order to determine if the phases can be successfully 
extracted with a small statistical error. Clearly, this simulation does not account for the details 
concerning the detectors. The possible impact of systematic errors will be discussed below. 

We have generated 10 5 events using a decay amplitude of the form in Eq. Q. As a first guess, 

(-) n 

we include in the decay only three channels, namely those mediated by the resonances D and 
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D*^, and the direct nonresonant decay — * (DfD°ir )j^ R . The form of the functions Aj for the 



Aj = BWAm 2 -) (p B ■ p n ) Fj (ml) , (10) 



resonances j = 1 , 2 is |Sj 

^3 ~ VY 3\""jJ \Ptl " Fit J ^jV'j 

where the invariant masses m| are defined as in Eq. (fljl. Fj(rrij) is a form factor, and BW 3 (s) is 
the Breit-Wigner function 



= — : rn > ( n ) 

(-) , 

mflj being the resonance masses = D , i?2 = -D^ )• For each j, the B and 7r meson three- 
momenta in Eq. (|1U|) are evaluated in the rest frame of the corresponding intermediate resonance. 

We have taken the usual expressions [3] for the form factors 3 and for the momentum-dependent 
width F R ,(s). The latter is given by 



(TIP. 



p(s) 



(12) 



where F Rj is the on-shell resonance width, and p(q 2 ) stands for the three-momentum of the res- 
onance decay products when the resonance mass is sfq* . The shape of the nonresonant decay 
amplitude, which is in general unknown ^3], has been taken — as it is usually done — as a constant 
function. In any case, as it is discussed below, this assumption has a negligible impact on our 
results. 

In order to carry out the generation of events, we need to introduce as input data the values 
for the physical quantities ay, 9j and the resonance widths. The expected relative weights aj of 
the two resonant channels can be obtained from the known branching ratios BR(B~ — ► D*°D~), 
BR(D*° -» L> vr°), BR(B~ -» D*~D°) and BR(D*~ -» D'tt ). We have [H] 



«i / BR{B~ -» D*°Dg) x BR(D*° -> D°7r°) ( 
~a~2~ \j BR(B~ -> D*~D°) x BR(D*~ -> Dsiv ) ' 

On the other hand, the nonresonant decay amplitude is uncertain; it is just expected to be smaller 
than the resonant channels. Taking into account that only the relative values between the three 
coefficients aj have a physical meaning in the simulation, we have taken a\ = 1 and «2 = 0.25, 
while for aN R we have considered different values, ranging from to 0.1. Concerning the phases 
0j, one expects CP-conserving parts to be relatively small, whereas the CP-violating SM phase 
is essentially the same for all amplitudes and can be factorized out. Thus, assuming that SM 
contributions dominate, it is reasonable to choose all phases 9j in our numerical simulation to be 
small numbers. In fact, it will be seen that this assumption is not relevant to our conclusions. 

Finally, other relevant inputs in our simulation are the on-shell widths T R of both resonances: 
events coming from a narrow resonance should be concentrated in a given region of the plot, 

hence the interference region between both resonances is expected to be relatively small. Present 

(-) , 

measurements of D* u and D* widths are not conclusive, giving in both cases only upper bounds 
of about 2 MeV. For our simulations, we have chosen to consider values ranging from 0.01 to 1 



3 The actual shape of the form factors in B decays is in general unknown. We have considered expressions similar 
to those used for D decays, finding that their incidence is not relevant to the discussion in this work. 
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MeV. The recent measurement of the D* + width, which is found to be around 0.1 MeV |14j . can 
be thought as a hint of the expected orders of magnitude. 



Figure 1: Dalitz plot for the B —* D s D°tt° decay 

As an example, we show in Fig. 1 the Dalitz plot generated with a^R = 0.1, 9\ = 0, 62 = 20°, 
&nr = 10°, and equal widths of 1 MeV for both resonances (besides the already given values of 
ai = 1 and «2 = 0.25). One observes that, even if both resonances are quite narrow, the events 
appear to be spread out in a large region of the plot. This is the consequence of a purely kinematic 
effect, due to the fact that both resonances are located very close to the threshold of the phase 
space. This effect compensates the narrow width suppression, and brings a good hope to extract 
the relative phases successfully. 

After carrying out this simulation of the decay, we have performed a fit of the data according to 
the fitting function given in Eq. where now the coefficients ctj and 9j are left as free parameters. 
In fact, as explained above, the fit provides only relative values for both amplitudes and phases [3], 
therefore we have kept fixed the reference values ol\ = 1 and 9\ = 0. The result of the fit is given 
in Tabled The method allows to extract the phase 62 with a statistical error as small as 1.4°. 



channel 


atj 


0.i 


D* U D S 


fixed 


fixed 


D°D* 


0.2514 ±0.0017 


(20.7 ±1.4)° 


nonresonant 


0.1007 ±0.0020 


(9.1 ± 1.2)° 



Table 1: Fitting results of the Monte Carlo sample. The events have been generated with 02 = 0.25, 
9 2 = 20°, a NR = 0.1, 9 NR = 10°, and r D „o = T D * = 1 MeV. 

We have performed a systematic study of the results of the fit allowing reasonable ranges 
of variation for the unknown quantities used to generate the Monte Carlo sample, namely the 
resonance widths, the weight olnr and the relative phases 0j. As a first outcome of this analysis, 
it is found that the statistical errors are independent of the initial values of the phases. Secondly, 
the errors for both the extracted amplitude and phase of the D* mediated decay (channel 2) are 
independent of the weight qnr of the nonresonant channel, even in the limit a^R = 0. This shows 
that the interference between the two resonant channels is not mediated by the nonresonant one, 
but arises from the above mentioned spread out of the events corresponding to resonance-mediated 
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decays. Finally, as expected, it is found that the errors in the extracted weights cej are independent 
of the resonance widths; on the contrary, the values of the widths do affect the quantity we are 
interested in, i.e. the error in the extracted relative phase 62 — 0\. This dependence is illustrated by 
the results in Table |2J where we have considered several simulations in which the amplitudes and 
phases atj, Oj have been taken as in the previously described example. We quote in the Table the 
errors obtained in the extraction of #2 — Q\ for different values of T D *o and Tp* . In the first five 
rows of the Table we have assumed equal D*° and D* widths, while in the last row we have taken 
T D *o = 100 KeV, ]?£>* = 10 KeV (in fact, a relative suppression of the D* width could be expected 
since the strong decay D* — * D^ir violates isospin). We see here that for a width as narrow as 
10 KeV the phase difference can still be extracted with relatively low statistical error. 



T D . ;T D * (MeV) 


Error 


1 


1.4° 


0.5 


1.5° 


0.1 


1.7° 


0.05 


2.3° 


0.02 


5.1° 


0.1 ; 0.01 


4.0° 



Table 2: Errors in the extracted value of 62 — 0i, for different values of resonance widths. Input 
amplitudes and phases for the event generation are same as in Table I. 

Before ending this section let us say a few words about systematic (experimental) errors in the 
extraction of phases. The evaluation of these errors is in general a quite difficult task. In order to 
carry out the complete analysis, one should perform a full numerical simulation of the experiment 
including the detector, which is out of the scope of this paper. Nevertheless, in order to have an 
estimate we can take into account the results from recent DP analyses jlj |5| . The latter suggest that 
the systematic error in the measurement of phases for intermediate channels with large branching 
fractions should not be above a few degrees, i.e. of the same order of those quoted in Tabled 

IV Expected size of new physics effects 

Let us now turn back to Eq. Q and analyze the theoretical expectations for the size of Ad^ m 
the context of a theory beyond the SM, in order to evaluate if this observable has potential chances 
to provide experimental evidences of new physics. To carry out the theoretical analysis we take 
into account the low-energy effective Hamiltonian relevant for the processes under consideration, 
including QCD corrections at the leading order. Then, to deal with long-range matrix elements, we 
use the simple factorization approach which should be adequate to estimate the significance 
of the new contributions |16j . 

In view of the large hadronic uncertainties and the usual amount of freedom to fix new physics 
parameters, we do not intend to perform an accurate calculation of possible nonstandard contri- 
butions to the — ► D^D ir° decay amplitude. Just as an illustrative example, we consider the 
rather representative framework of multihiggs models, showing that the situation becomes quite 
promising if nonstandard contributions to penguin diagrams are comparable to those arising from 



8 



SM physics. 

Our theoretical analysis is based on the AB = 1 effective Hamiltonian |1U1 117j 



n cS = -j= \ v cb vZ(c 1 o 1 + c 2 o 2 ) - v tb v t * s (J2 CM) 




8=3 



(14) 



where C, are Wilson coefficients evaluated at a renormalization scale \i ~ m^, and Oj are local 
operators, 



Here V ± A refers to the Lorentz structure 7^(1 ± 75), a and /? stand for SU(3) color indices, T?g 
are generators of SU(3) color transformations and G a ^ v denotes the gluonic field strength tensor. 
Contributions from electroweak penguins will not be taken into account, therefore these operators 
have not been included in (|14|). We will also neglect the effect of the electromagnetic dipole 
operator. Within the SM, the coefficients Cj can be calculated at the scale mw, and then evolved 
to fi « nit, through the renormalization group equations ^Uj. The Vckm factors corresponding 
to each operator have been explicitly separated in (|14[). so that with good approximation the 
coefficients in the SM can be assumed to be real numbers 4 . Moreover, in view of the unitarity 
of the Vqkm matrix, one has V^V^* = —V c bV* s — V u bV* s ~ —V c bV* s , where the correction due 
to the V u bV* s term is about 5%. In this way, for the case under consideration, the CP-violating 
phase carried by the penguin contributions in the SM is approximately the same as that coming 
from the tree operators C\ and C2, and will factorize out for the decay amplitudes of interest 
(the contribution of the V u bV* s term to the full amplitude will be below 0.5% if, as expected, the 
total penguin amplitude does not exceed 10% of the tree piece). In a given extension of the SM, 
however, the coefficients C{ will carry in general nonvanishing CP-violating phases ipi, allowing for 
the interference effects discussed in the previous sections. 

In general, in a theory including physics beyond the SM, one expects that the new particles 
can be integrated out at the mw scale, leading to new contributions to the coefficients Cj(mjy)- 
However, since the new particles have been integrated out, the running of the coefficients down to 
fi ~ nib proceeds just as in the SM ^Sj. This running of SM coefficients has been analyzed in detail 
in Refs. I19j and will not be repeated here. 

In the evaluation of the amplitudes (VP\TC c f[\B) , the scale and renormalization scheme depen- 
dence introduced by the coefficients C% should be compensated by that of the matrix elements of 
the quark operators Oi between the hadronic states. However, as stated above, to evaluate these 
quantities we will use the factorization ansatz, and in this approach the matrix elements are written 
in terms of decay constants and form factors, which are both scale and renormalization scheme in- 
dependent. In order to achieve the required cancellation, it is possible 47^ to calculate the one-loop 
corrections to the partonic matrix elements (scc\Oi\b), and to define new effective coefficients Cf 

4 In fact, they carry small CP-violating and CP-conserving phases, coming from Cabibbo-suppressed contributions 
and absorptive parts of loop diagrams respectively. 



Oi 

o 7 



(c Q b a )v-A(s/3 Cp)v-A 
(s a b a ) V ^ A Eg' % QpW+A 

{g s /%Tt 2 )m b s a a** (1 + 75) T%* bp G 



'a 



O2 
6 



(c/3 b a )v-A(s a Cp)y-A 

(sp b a)v-AY, q >{Q' a q'f3)v-A 

(spbaW-AY.q'iQaq'pW+A 



(15) 
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such that the one-loop quark-level matrix elements read 



(scc\n e{i \b) = Y,Cf(scc\O l \b) tree . (16) 



At NLO these effective coefficients will be given by the original plus QCD corrections, 

cf = cm + g J2 K MCM . (17) 

The analytic expressions for the functions Kij can be found in Refs. |17l 1201 |2"T] . It can be shown 
that now the effective coefficients Cf are scale and scheme independent, as well as gauge invariant 
and infrared safe [22]. An important point is that the corrections introduced in Eq. (|17|) involve 
the coefficient C7, which can receive important contributions coming from nonstandard physics, as 
occurs e.g. in the case of two-Higgs-doublet models ^Hj. Even if the operator O7 does not contribute 
directly to the B — ► VP decay amplitudes in the factorization approach, the combination in l|17jl 
implies that the new physics corrections to C-j are translated to other effective coefficients Cf s and 
thus to the decay amplitude. 

The previous analysis can be now applied to the decays of our interest, namely the resonant 
processes B~ -► D*~D°; D*~ -> Z^vr and B~ -► D*°D-; D*° -> D°tt° that dominate the 
three body decay B~ — > D°D~ir°. In the described framework, the relevant two-body amplitudes 
(D*- D°\H cS \B-) and (Z>*°.D7|?4 ff |p-) will be given by 

6 

(VP\H eS \B-) = -^V^V^CfiVPlOilB-}^ 

i=l 



^=V cb V; s a(B~ - VP)X ^-^ p ) , (18) 



where the subindex FA denotes that the matrix element is evaluated within the factorization 
approximation. The factor a(B~ — ► VP) includes the effective coefficients Cf s , whereas X^ B ~^ VP ) 
contains the form factors related to the factorized amplitudes. For the processes under consideration 
one has [231 HI 



2 

TTln 

a{B~ -> D D~) = a x + a 4 - 2a 6 



(m& + m c )(m s + m c ) 
a(B- -» D*~Z) ) = ai + a 4 , (19) 

where the coefficients a* are defined as a, = Cf+CffjiNf),, for i = 1, and a, = Cf+Cf 1 /(Nf) i 
for i = 4,6. The effective parameters (N° s )i in these expressions account for the uncertainties 
introduced when calculating the matrix elements of the effective operators between hadron states 
dll HOI HI]. The factors X^ B ~^ VP ^ are given by 

x (b--.b*o D7) = 2 f Ds m D *oA BD \ml s ){e* D * .P B ) 

x (B->nr D °) = 2/n. m n.Ff D (mL)(4.P fl ), (20) 
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where Ey are the corresponding V meson polarizations, Pg is the B four-momentum, and the 
expressions include decay constants and form factors that can be estimated in specific models. In 
fact, Eqs. (|20j) have been quoted only for completeness, since the factors X^ B ~^ VP ) cancel out 
in our estimation for A0^ exp \ This can be seen by noticing that the expression for Ad^ in 
Q involves ratios between SM and new physics amplitudes. According to previous assumptions, 
the effects of new physics are only present in the effective coefficients Cf — or, equivalently, 
a{B — > VP) — , therefore any global factor will cancel. One has in this way 

a NP e -iv NP (D*°D-\H cS \B-) NP ( ai + a A -2ra 6 ) NP 



a SM (D*QD7\H eS \B-) SM (a 1 +a 4 -2ra 6 ) SM 

^ p e~^ NP {D* s -D°\7i c{i \B~) NP {a x +a A ) NP 



al M {D* s ~ D0\H eS \B-)SM - ( ai + a4 )SM 



(21) 



where r stand for the mass ratio m 2 D /[(mb + m c )(m s + m c )], and — as in the previous sections — we 
have assigned labels 1 and 2 to the channels mediated by the resonances D*° and D* respectively. 
Average values of quark masses yield r ~ 0.5. 

In order to analyze the possible NP effects in our observable A9^ exp \ let us consider the typical 
situation of a theory including an extended scalar sector. In the case of multihiggs (MH) models, the 
scalar-mediated tree contributions to C\ and C2 can be neglected, since in general scalar couplings 
are proportional to the current quark masses of the involved vertices. On the other hand, penguin- 
like diagrams mediated by the new scalars involve vertices which are proportional to the top quark 
mass, thus they are potentially important. Then, while SM amplitudes are dominated by tree 
contributions (Cff > Cf M for % = 3. . . 6), in a MH scheme the main effect of the extended scalar 
sector on a\ and Q2 occurs through the new contributions to the effective coefficients 04 and ciq. 
In this way, from Eqs. © and (|21|) one gets 

/ n MH n MH\ \n MH \ 

A9^ * 2 sin^ ( V " %W J ~ "4r sin/^^J . (22) 

V Off «2 / 01 

As a first outcome from this expression, it is seen that the ratios a^ p /otj do not cancel with 
each other, consequently the asymmetry A6^ exp ^ is in principle nonzero. 

Even if the result in (|22|) is just an estimate, it can be taken into account in order to show 
that new physics effects can be significant enough to provide an observable signal. According to 
the analysis presented in the previous section, this would be achieved if new physics contributions 
to a6 reach about 10% of the SM tree amplitude, and carry a CP-violating phase of order one 
(one would obtain in this case an asymmetry A0( ex ri of about 10 de grees). Within the SM, the 
effective coefficients |ai| (tree) and \ae\ (penguin) are estimated to be approximately 1 and 0.06, 
respectively |23] . Thus, one would have important chances of measuring nonstandard physics if 
new contributions to a§ carrying large CP-violating phases are comparable in size to SM ones. It 
is worth to point out that this level of contribution of nonstandard physics is indeed suggested by 
some puzzling experimental results on penguin-dominated modes, such as the B — > r/'K branching 
ratios |25| and the time-dependent CP asymmetries in B° — > 4>K$ |26| . The experimental values 
for these observables are at least 2a away from SM expectations, and can be seen as indications of 
large new physics effects at the penguin level. 

We believe that these experimental observations on penguin-dominated B decay channels al- 
ready provide a substancial ground to encourage the DP analysis of B^ — > D^D ^ proposed here. 
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On the other hand, we point out that the room for nonstandard contributions to penguin ampli- 
tudes — and thus to the phase difference A9^ exp ^ — is relatively large, mainly due to the existing 
theoretical uncertainties in the evaluation of SM amplitudes, and to the large number of unknown 
parameters included in most scenarios of new physics. To be definite, let us take here as an exam- 
ple one of the simplest possible extensions of the SM, namely a two-Higgs-doublet model (THDM) 
type III. In particular, we consider a minimal scenario j^Tj which does not include tree level FCNC, 
and the number of new parameters is reduced to four (three Yukawa couplings parameters plus 
the charged Higgs mass). In this framework the main new contributions to b quark decays arise 
from one-loop diagrams involving a virtual top quark, while neutral Higgs-mediated diagrams are 
shown to be negligible |181 I28| . As stated, in this kind of models the largest new contributions 
to the amplitudes <Zj come through the dipole coefficient C7, and the allowed space for the new 
parameters is mainly constrained by the effects on B — > X s ~j decays |28| . Taking into account the 
bounds in Refs. [27|[2Hj, it is possible to estimate the allowed values for both the amplitude a§ 
and the CP-violating phase (p. We find that within this model the phase difference A8^ expS) can 
be as large as 3 degrees, which, according to the analysis Sect. Ill, would be around the limit of 
observability for the number of events considered. 

The example below should be taken just as an illustration to show the potentiality of our analysis 
through a simple manageable case. Clearly, the inclusion of more degrees of freedom would relax 
the experimental bounds on the new model parameters (imposed e.g. by the chosen mechanism 
to avoid unwanted flavor changing neutral transitions), allowing higher values for the measurable 
phase difference A0( exp ^ which will exceed the observability limits. In addition, other possible 
frameworks of nonstandard physics have been shown to provide enhancement effects on penguin- 
dominated processes, offering an explanation for the puzzling time-dependent CP asymmetries 
in B° — > cpKs- Among the most popular scenarios, recent analyses include R-parity violating 
supersymmetry [212], left-right supersymmetric models [201; an d theories including warped extra 
dimensions [3J. In all these models — which include in general a rather large number of new 
parameters — , it has been shown that new physics contributions can be of the same order as SM 
penguin amplitudes. In this way, their effects on the b — > ccs channel could provide an observable 

signal in the DP analysis of — > D^D°tt decays proposed here. 

V Summary 

We discuss the possible measurement of nonstandard CP violation in B^ — > D^D ir°, exploiting 
the fact that for these processes the asymmetry between B + and B~ decays is expected to be 
negligibly small in the Standard Model. The presence of two resonant channels provides the neces- 
sary interference to allow for CP asymmetries in the differential decay width, even in the limit of 
vanishing strong rescattering phases. 

In order to measure the CP-odd phases entering the interfering contributions to the total decay 
amplitude, we propose to use the Dalitz Plot fit technique. This allows a clean disentanglement 
of relative phases, independent of theoretical uncertainties arising from FSI effects. The expected 
quality of the experimental measurements has been estimated by means of a Monte Carlo simulation 
of the decays, from which we conclude that the phases can be extracted with a statistical error 
not larger than a couple of degrees, provided that the widths of the intermediate D*° and D* 
resonances are at least of the order of a hundred keV. On the theoretical side, within the framework 
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of generalized factorization we perform a rough estimation of possible nonstandard CP violation 
effects on the interfering amplitudes. We take as an example the typical case of a multihiggs model, 
showing that the level of accuracy of the DP fit measurements can be sufficient to reveal effects of 
new physics. 

Let us finally stress that tree-dominated decays like — * D^D°tt° are usually not regarded 
as good candidates to reveal new physics, since the effects on branching ratios are not expected to 
be strong enough to be separated from the theoretical errors. Our proposal represents a possible 
way of detecting these effects by means of CP asymmetries, which can allow the disentanglement 
of new physics contributions to penguin-like operators in a theoretically simple way. 
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